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Abstract 

A detailed description of the method for analytical evaluation of the three-loop 
contributions to renormalization group functions is presented. This method is employed 
to calculate the charge renormalization function and anomalous dimensions for non- 
Abelian gauge theories with fermions in the three-loop approximation. A three-loop 
expression for the effective charge of QCD is given. Charge renormalization effects in 
the S'[/(4)-supersymmetric gauge model is shown to vanish at this level. A complete 
list of required formulas is given in Appendix. 

The above-mentioned results of three-loop calculations have been published by the 
present authors (with A.Yu. Zharkov and L.V. Avdeev) in 1980 in Physics Letters B. 
The present text, which treats the subject in more details and contains a lot of cal- 
culational techniques, has also been published in 1980 as the JINR Communication 
E2-80-483. 



1. Introduction 



The renormalization group method when apphed to asymptotically free models 
results in an "improved" perturbation theory. Its expansion parameter, an effective 
charge g'^{Q'^/A'^,g'^), decreases logarithmically with the increase in the momentum 
transfer Q^. The existent QCD calculations of various deep inelastic processes in the 
first two orders in appear to be consistent with the present experimental data [1]. 
However, the next-to-leading corrections (i.e., those ~ ) are fairly large. It leaves 
open the possibility that the higher-order contributions will be important. 

The calculations in higher orders are also of interest from another standpoint. They 
might serve us a starting point for summing the perturbation theory expansions of 
QCD, as it is done, for instance, in the (p'^ model [2]. Moreover, these calculations 
can shed light on some peculiar aspects of certain field theory models. For example, in 
the 5'C/(4)-supersymmetric non-Abehan gauge model derived in [3,4], the charge renor- 
malization effects are shown to vanish to the two- loop order [5]. The corresponding 
three-loop calculations presented below give the same answer: The charge renormaliza- 
tion function (3{g'^) is equal to zero. Apparently, the vanishing of (3{g'^) at the three- loop 
level is not a sheer coincidence, but an indication that this effect holds to all orders. 

The first three- loop QCD calculation in the framework of the renormalization group 
has been performed in [6], where the total cross section of the e~^e~ - annihilation into 
hadrons has been computed analytically. This result is confirmed in [7] by a numerical 
calculation and in [8] also analytically. However, these calculations involve the (3{g^) 
function to order g^, whereas all other three-loop QCD calculations require the next, 
~ g^, contribution to /3{g'^). The charge renormahzation function /3{g'^) for the non- 
Abelian gauge theory including fermions is known to g^ only, i.e., in the two-loop 
approximation [9]. In the present paper we describe a method which enables one to 
evaluate (3{g^) at the three-loop level. We present the results of these calculations and 
the full list of needed formulas. 

2. Renormalization group in the minimal subtraction scheme 

We consider a non-Abelian gauge group theory with fermions belonging to the 
representation R of the gauge group G: 

11 ^ - 

m=l 

Gl, = d,Al - d^Al + gr'-^A^Al , P.^f = d.i^T - wRtji^TK ■ 

Here 77" is the ghost field, a is the gauge parameter, and /"^"^ are the totally anti- 
symmetric structure constants of the gauge group G. The indices of the fermion field 
tp"^ specify color (i) and flavor (m), respectively. The matrices i?" obey the following 
relations: 

[R'',R^\-^ir^R'', fo^fbcd^ fj^gab^ R^R^^GfI, tr(i?"i?^) = (2) 



2 



In particular, the values of group invariants Ca, Cp and T in the fundamental (quark) 
representation of SU{N) are: 



2N ' 



T 



1 
2 



(3) 



The underlying gauge symmetry of the Lagrangian ([T]) gives rise to the well-known 
Slavnov- Taylor identities [10] extensively used throughout the paper. In particular, 
a transversality of the radiative corrections to the gluon propagator allows one to 
compute such a correction in the scalar form, i.e., with its Lorentz indices contracted. 

We now turn to a brief discussion of the renormalization procedure. In this paper 
we adopt the renormalization prescription by 't Hooft [11], the so-called "minimal 
subtraction scheme", which by definition subtracts only pole parts in e from a given 
diagram. The renormalization constants relating the dimensionally regularized IPI 
Green function with the renormalized one, 



R 



/i2 



lim Zi 

e-i>0 



] r {Q^,aB,gl,€) 



look in this scheme like 



Zr(l,a,gA=l + f2cPic^,9')e 

^ ' n=l 



(4) 



(5) 



with e = d being the space-time dimension. In /i is the renormalization 
parameter with the dimension of mass. The bare charge g]^ is to be constructed from 
appropriate Z's. The most convenient choice is as follows: 



gl = ii'^g^~ZlZ^'~Zf. 



(6) 



Here Zi is the renormalization constant of the ghost-ghost-gluon vertex, Z3 and Z^ 
being those of inverted gluon and ghost propagators, respectively. Note also in (jH) to 

be given by = oiZ^. The Green function F/j a, g^^ satisfies the renormalization 
group equation 



(9 (9 



dg' 



d 

j3{a,g^)a— - 7r(a,5'^) 



R 



-^,a,g 



(7) 



and the normalization condition F/j (^^,a,0^ = 1. The anomalous dimensions 7r are 
given by the relation 



,2 d ji) 



dg 



-c^'[a,g 



Similarly, from 



9l 



2e 



n=l 



(9) 



one obtains the charge renormalization function /3 
> d 



Kg' 



9 



dg' 



1 a^'Kg') = 9' [27i(«,/) - 73(«,/) - 273(a,^?')] , (10) 
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which is known to be gauge independent [12]. Thus, the computation of 7r(a,5'^) 
and P{g^) requires the functions c^j^\a,g'^) for the renormahzation constants in the 
right-hand side of (jH]). 

The residues of higher-order poles in the expansion ([S]) and are related with c^^^ 
and a^^^ by the equahties 



d d 



9 



cy{a,g 
d 



2 9 {n+l), 2\ 



:iii 



(12) 



We choose to work in the Feynman gauge a = 1 throughout this paper. For checking 
the higher residues by means of flTT]) one may use the results of the corresponding 
two-loop calculations [13] performed in a general gauge. 

According to the minimal subtraction prescription [11], the renormahzation con- 
stants are uniquely determined by requiring that all the divergences in e disappear 
from the product Zy 5'^) T (Q^, a^, e), so that the limit e — )■ in does ex- 
ist. However, we find a somewhat different (but equivalent) definition [14] to be more 
convenient: 

Zt = 1- ICR'T. (13) 
An operator /C picks out all the pole terms in e, 



(14) 



n<0 



R' is the BPHZ minimal subtraction procedure (/^-operation) with its final subtraction 
missing: R = (1 — /C)i?'. In other words, the /^'-operation subtracts all the divergences 
of internal subgraphs but does not subtract an overall divergence of a diagram. To 
construct R' explicitly one can employ the following recursion relation [15]: 



B!G = G + Y^ i-JCR'Gi] 



-]CR' Gm) ■ G /{Gi -\- ... -\- Gr, 



(15) 



where the sum is over all sets of disjoint IP/ divergent subgraphs of the diagram G, 
and G/{Gi + ... + Gm) is the diagram obtained from G by contracting Gi, ...,Gm to 
points (as an example see Fig.l). 



R' 





KR''^' 




-2 



KR' 





R'' 




KR' 




Fig. 1 
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The KR'G is the negative of a contribution from G to an appropriate renormahza- 
tion constant. The computation of JCR'G is simphfied drastically owing to the following 
fact [16]. 

Let a diagram G be infrared finite in a range of external momenta ki and internal 
masses rrij. Then in this range ICR'G is a polynomial in ki and rrij. Therefore, it either 
is independent of ki and mj (for a logarithmically divergent diagram G) or loses such 
a dependence after differentiating once or twice with respect to ki. 



3. A method for computing three-loop integrals 



This feature of JCR'G provides the basis for a simple and efficient computational 
technique developed in [15], which enables one to evaluate analytically all three- loop 
contributions to the renormalization group functions 7 and (3 in any renormalizable 
theory. It is shown in [15] that one may calculate ICR'G (properly differentiated, 
if necessary) with all its external momenta equal to zero and with an auxiliary mass 
m 7^ introduced into one of its internal lines (which is sufficient to remove all infrared 
divergences). The momentum integration corresponding to this line is chosen to be the 
last one. It looks like 

* (16) 



/ 



(p2)"(p2 _^ 

and is readily done using Eq. (166|) in Appendix. We thus show the last momentum 
integration to be trivial. Therefore, the problem of three-loop calculations reduces to 
computing the two-loop massless integrals depending on a single momentum p"^, 



dt dq 



t^<^q^P{p-tf^{p-qY'^{t-q) 



2p 



(17) 



with a, 13 
p = 0,-l. 



7, a and p being integers. If one of the denominators is missing (e.g., 
■2, ... ) the integral (|T7|) can be evaluated by sequential use of Eq. (EZj). Oth- 
erwise one needs the non-trivial two-loop integration formulas deduced in [17] through 
the x-space Gegenbauer polynomial technique. In Appendix we give a list of relevant 
integrals of the type f lT7|) . 

As an illustrative example we consider an integral 



J 



dp dq dt (qt)' 



p^qH'^{p — qY{p — ty{k — qy{k — ty 



:i8) 




Due to quadratic divergence, it should be differentiated twice with respect to k. 
Using the relation 



92 



dkf, dkf. 



{k-qY{k-ty 



8{k -q){k-t)+ 4e[{k - qf + (A; - tf] 
{k-q)\k-tY 



(19) 



we obtain JCd^R'J as displayed in Fig.2 in self-evident notation. Since ICR' J = k'^A (i) , 
we finally get 



All 



1 



8-Ae 



-Kd'R'J 



[in 



2\3 



24^2 



1 

32e 



(20) 



The last two diagrams in Fig.2 diverge logarithmically so that one can compute them 
with = provided that a non-zero mass is introduced into one of the differentiated 



lines, i.e., into that with a blob. 






R' 




-2 




-2 



KR' 





Fig.2 

The problem of evaluating KR'G at the three-loop level thus reduces to the integrations 
f lT^ and f lT7|) . The described procedure has been employed in a considerable part of 
the calculations presented in this paper. 

One can also determine the pole part of (|T8|) . /CJ, by means of a somewhat different 
method, which involves transferring an external momentum to the other vertex in order 
to simplify the denominator. 





Fig.3 



Consider the difference (Fig.3) 

dp dq dt {qtY 



I 



1 - 



{k-q) 



21 



p^qH^ip - qY{p - tY{k - qY{k - tf 
dp dq dt {qtY{2k^q^ — k"^) 

pH^q\p - qy{p - ty{k - qy{k - ty 



= '^k^J^ — k'^Ji . 



(21) 



Let us further subtract from the other integral having a more simple structure of 
the denominator: 



J„- 



dp dq dt q^{qty 



dp dq dt q^{qty[2k^ty — k"^] 
'^^ ' pH^q^p - qY{p - tY{k - qf ~J p^qH^p - qf{p - tf{k - q)^{k - tf 
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• (22) 



There is only one (logarithmically) divergent integral in the right-hand side of fl22|) . 
namely 

dp dq dt 2q^ty{qtY ^^^^ 



I 



p'^qH\p - qY{p - tY{k - qY{k - tf ' 

Due to the absence of divergent subgraphs, its pole part does not depend on k and 
coincides with 

/■ dp dq dt 2q^t^ {qtf 

As to the integral Ji, it diverges logarithmically and contains divergent subgraphs. We 
note the difference 

J _ f t/p dq dt {qtf 

' J p^qH\p - qfip - tY{k - qf ^ ' 

to be convergent, and combining the last five relations finally obtain 

^j_r f dpdqdt {qtY[A{kt)ikq) + 2qH' - t\k - qf] _ 1 25 

J p^qH^ik - qY{p - qY{p - tf ^'"^ ' \ 128^ 32£ 

(26) 

This integral is easy to evaluate with the use of formulas listed in Appendix. Adding 
to fl26p the appropriate counterterms gives for ICR' J the same answer as in (I2U]) . 

The essence of the procedure presented above is as follows. One subtracts from the 
initial integral J an infrared finite integral J' with a more simple denominator reducing 
thus the degree of divergence. Such a subtraction is to be repeated until the difference 
becomes convergent. 



4. Calculation of specific diagrams 

It is now seen that the three-loop momentum integrals contributing to Z's are al- 
ways calculable. However, one must introduce an auxiliary mass into the diagram 
(which as a rule represents a sum of distinct integrals similar to (|T8|) ) and into all its 
counterterms in a consistent fashion. For the most complicated diagrams of the gluon 
propagator this task appears to be unmanageable. Therefore, we deal with the dia- 
grams of the topological type, depicted in Fig. 4, as follows. We reduce the numerator 
of the integrand to the scalar form and then decompose it into a sum of invariants like 
k'^{q — p^q^ip — t)"^, .... Canceling numerator against denominator and taking sym- 
metry into account results in at most 66 distinct three-loop massless integrals. Their 
pole parts are to be found either by direct use of fl671) - fl72l) or by differentiating, intro- 
ducing a mass, and then converting ICR' into /C through the compensating subtraction. 
The latter pole parts are given in Appendix. 




Fig.4 
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# <5> €> €> O 



Fig.5 

The propagator diagrams of more simple ("nested") topology (Fig.5) can be com- 
puted straightforwardly using ( I67|) - (l72l) . The remaining topological type is represented 
by a single diagram (all others equal zero owing to the antisymmetry of the group 
structure constants) which can be easily calculated by means of differentiation: 



,e^(C,-C.)(C.-^)/16^20_32^pj^^^ (27) 




All the diagrams of the ghost-ghost-gluon vertex diverge logarithmically. We eval- 
uate them setting all external momenta to be zero and introducing an auxiliary mass 
into one of the internal lines. For each particular diagram this "potentially infrared" 
line is easy to identify. 

Thus, all the diagrams of a certain Green function are calculated in the same fashion: 
with an auxiliary mass for the vertices and without it for propagators. It enables one 
to perform the subtractions either following 't Hooft [11] or determining JCR'G for each 
individual diagram. In order to check the intermediate results we choose the latter way. 

The problem of evaluating the group weights appear to be of no substantial diffi- 
culty. Mostly it reduces to making contractions in the products of several structure 
constants f"'^'^. The following graphical representation is here of great use [18]. 



/ 



abc 
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[-Ca) => f^^^fb^ = -CaS'''' 

= > f'i--f-bgpcd^_l(j^fabc ^28) 
jgai jijd jjbh jheg jdce _ g _ 

The last two relations are derived from the Jacobi identity 

I ^^'^^^ ^ Jabcjade _j_ Jabejacd _j_ Jabdjaec g (29) 




The only products of structure constants which cannot be contracted by the sequential 
use of fl25]) are the following (Fig. 6). 





Fig. 6 



From fl2^ we obtain 





1 r'S 

'8 




(30) 



However, one fails to express the graphs of Fig. 6 separately in terms of Ca- In a specific 
case of the SU{N) group, we have found 




2 



(31) 



Fortunately, the relation ( l30i) is quite sufficient for the three-loop calculations of the 
renormalization group functions. Only the sum of the diagrams of Fig. 6 contributes 
to the final answer. This fact is easy to explain. The non-trivial products (Fig.6) 
might contribute to the vertex anomalous dimension, 71 (a, ^f^), only. But it is known 
to vanish in the Landau gauge: 7i(0, (7^) = 0. Hence these products do not contribute 
to the gauge independent function /3{g'^) and consequently, to 71 (a, 5^^) in arbitrary 
gauge as well. 

Concluding this section we wish to discuss one more example where Slavnov- Taylor 
identities [10] have been fruitfully used. To facilitate the computation of the vertex 
diagram with the two-loop three-gluon insertion 



JC 



1 



dp 



4:-2eJ (27r) V(P^ + ^^) 
we employ an identity 

P,rf:,{k,q,p) = G{p 

where a notation is as follows: 



R 



Mf;{k,q,p)V-\q'){q^g^, - q^q,) + \ p ^ 




(32) 



(33) 
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q.Mf;{k,q,p), 



P 



fx 



P 



V 



Sab 



(9,. - ^Mp') + a 



p^ 



P^lPv 

p2 



(34) 



In our case A; = so that ( l33|) transforms into 

P,T;11{0, -p,p) = G{p')V-\p'){p'g,, - p,p^)Mf;{0, -p,p) . 



(35) 



Identity (l35l) allows us to calculate M^^ rather than fairly complicated three-gluon 



vertex r^[^^ . 



5. Three-loop results for QCD 

A total number of topologically distinct three-loop diagrams contributing to f3{g'^) 
amounts to 440 (without counting opposite directions of the ghost and fermion lines). 
For performing the Lorentz and Dirac algebra, reducing the integrands, decomposing 
the scalar products, evaluating and summing standard integrals, the computer program 
SCHOONSCHIP [19] has been substantially used. The total execution time is rather 
difficult to estimate. Here we only indicate that the diagrams of Fig.7 require 110 and 
90 seconds, respectively, at the CDC-6500 computer. 





Fig. 7 
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Our final results obtained in collaboration with A.Yu. Zharkov are as folows (/ is 
the number of flavors, h = 



7i(l, h) = -^h - + (-^Cl + ^ClT/ ) , 

73(1, h)=h ih^A - \Tj\ + Ir? {^^C\ - 5CATf - ACpTf 



/4051 
\1U 



- oC(3) C 



875 

Is 



+ 18C(3) CiTf 



- ( ^ + 24C(3)^ CaCfT/ + 2ClTf + ^C^TV + "^CpT^f 



(36) 



(37) 



^ + 9C(3)) ClTf + + 12C(3)) CAC^Tf - ^CaT'/' 



,(38) 



+ '^CAC^Tf - fCpT'f - 2ClTf\ . 



(39) 



The cancellation of the transcendental ^(3) in the expression for /3(/i) is in complete 
analogy with QED treated in the minimal subtraction scheme, where [20] 



^Qi.D(a) = -- + 4^--^. 



(40) 



In a particular case of QCD, when fermions transform according to the fundamental 
representation of S'f/(3), [3{h) reads: 

^ <-xi . ?/) .<-xo. . |/) .;.(-^ . 555, _ . 

Now we are in a position to find an effective charge h {^,h^ from 



ln^= r ^ = ^(h)-m, 



where ijj(h) represents an indefinite integral j'* Let us express h in terms of 

renormalization group invariant quantity In ^ + ■0(/i) = In ^ = L, where A is the 
momentum scale. Assuming 



/^^ Jh 



(42) 



(43) 



11 



we arrive at 



PI 



(44) 



and obtain from (I32 
h{L) 



1 /3ilnL 5/3^ - /3i In /3o PlXn^ L InL 



+ 



1 



L3 



[/32/3o - Pi + /3i(5/3o - Pi In /3o) + (5/3o - Pi ^^Pof] + O 



In^L 



(45) 



with 5 being an arbitrary constant. Fixing the momentum scale A by choosing, as 
1 Pi\nL Pl{\n^ L-\nL) P2P0 - Pf 



usual, 5 = ^i^, we finally get 

Po 



KL) 



PoL PI L2 



+ 



PIL' 



+ 



PlL' 



+ 



In^L 



(46) 



Using fHTl) . fj43|) and fH6|) one readily finds the QCD effective charge in the three- loop 
approximation. 



6. Vanishing of (5{g^) in a supersymmetric gauge model 

Some time ago a very interesting SU (4)-supersymmetric non-Abelian gauge model 
has been derived [3,4] which exhibits the vanishing charge renormalization effects, since 
its charge renormalization function P{g^) proves to be zero through the two-loop order 
[5]. The Lagrangian is [4]: 




with a,b,c = 1, A^^ — 1; m, n = 1, 4; r,t = 1, 2, 3. Here £ym is the pure Yang- 
Mills Lagrangian with SU (N) gauge symmetry. The matter fields (Majorana spinors 
A^, scalars 0" and pseudoscalars Xr) transform according to the adjoint (regular) 
representation of ^[/(A^). Hence 

with similar expressions for and 'D^Xr- The six real antisymmetric 4x4 matrices 
a^, P"^ obey the relations 

= [P^,P% = -W\ \a\P'\. = 0. (48) 

The other properties of these matrices and their explicit form are given in Appendix. 

To determine the contributions to the renormalization group functions of the model 
( H7|) from the diagrams without scalar and pseudoscalar particles, one may use the 
results (I36])-(I39]) with 

Ca = Cf = N, Tf = 2N. (49) 
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This leads to 

/3(/i)without scalar^ = 'Nh^ + lOiV^/l^ + ^N^h\ (50) 

Now an appropriate scalar contribution must be added to (l50l) . In the two- loop ap- 
proximation it has been done in [5] with the intriguing result = 0. 

The method of our three-loop calculations is described above. Here we shall only 
consider the issue of applicability of the standard dimensional regularization to super- 
symmetric theories. This subject has been discussed by various authors [21]. Proceed- 
ing in the spirit of Ref. [21] we write down the following rules of the "supersymmetric 
dimensional regularization" which is to maintain both gauge invariance and global su- 
persymmetry: The relations defining the Dirac matrices look as in four dimensions (see 
Appendix) while the numbers of scalar and pseudoscalar fields equal 3 + e rather than 
3. This modification of the regularization maintains equal (and integral) total numbers 
of Bose and Fermi degrees of freedom even in 4 — 25 dimensions: 8 components of four 
Majorana spinors correspond to (2 — 2e) massless vectors -|- (3 + e) scalars + (3 + e) 
pseudoscalars = 8 bosons. It is this matching of the Fermi and Bose field components 
that is crucial for preserving supersymmetry [21]. 

For lack of a rigorous proof, we have verified the invariance of the supersymmetric 
dimensional regularization by direct calculation of (5{h) at the two-loop level in two 
different ways: 

m = mAh)-i,{h)-2'^^{h)] (51) 

and 

P{h) = /i[274(/i) - i^{h) - 27a(/i)]. (52) 

Here 71 and 74 are the anomalous dimensions of the ghost-ghost-gluon and fermion- 
fermion- scalar vertices, and 73, 73, 7,/, and 7a are those of gluon, ghost, scalar and 
fermion propagators, respectively. In the standard (with 5^'' = 3) dimensional regular- 
ization, these anomalous dimensions are (in the Feynman gauge): 

11 = -^- lN'h\ 74 = -5Nh + 5N'h\ 

^3 = -2Nh + = -2Nh, (53) 

73 = ^ - N^h\ 7a = -4:Nh + 6N'^h\ 

With the use of supersymmetric dimensional regularization (with S^"^ = 3 + e), we 
obtain 

71 = - lN^h\ 74 = -5Nh + ^iV^/^^ 

73 = -2Nh + N'^h^, -f^ = -2Nh - N^h^, (54) 

73 = ^ - iN'h', 7A = -^Nh + QN'h'. 

Using flSTl) gives P{h) = for both regularizations while fl52|) leads to P{h) = —2N^h^ 
for the standard regularization and to P{h) = for the supersymmetric one. This 
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discrepancy shows the former regularization to be noninvariant under supersymmetric 
transformations. 

For our three- loop calculations we employ formula fl3T]) . Below we write down the 
scalar contributions to anomalous dimensions through the three-loop order calculated 
in the supersymmetric dimensional regularization scheme (in collaboration with L.V. 
Avdeev) : 

= -Nh + ^N^h^ + (f - i C(3)) N'h\ 
7r' = - fN^h^ + + l C(3)) N'h\ (55) 

-seal ^ ^N^h\ 



From (1551) and flSTl) we obtain 



32 



/3'^''\h) = Nh'^ - lON^h^ - ^N^h^ (56) 
and using (150|) . arrive at the final result 

/3{h)three loops = 0. (57) 

It is worth mentioning that the use of the standard dimensional regularization yields 
7^ = -Nh + ^-^N^h' + - i C(3)) N^h\ 

ir' = -Y^'h' + + l ^(3)) ^r' = i^'^'' (58) 

(3 (h) three loops = SN^h'^. 

The result (157|) implies the absence of the charge renormalization effects in the 
model (H7j) to the three- loop order. It confirms a conjecture that (3{h) in this model 
vanishes to all orders. If it were the case, the model fj47|) would be unique in the four 
dimensional quantum field theory. The vanishing P{h) might imply, for instance, that 
this model would be free of supersymmetric anomalies [22]. In any case, a rigorous 
argument proving this conjecture on symmetry ground is now a great urgency. 

We would like to thank L.V. Avdeev, G.A. Chochia and A.Yu. Zharkov for the help 
in some calculations. 

APPENDIX 
Feynman rules for the model (1) 
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T] --* Tj 



- — 8°-^ 



P 



2 V ' 



P 



C 



P^ 



n 




gp.r'" 



rrv I 




^P 



/3 /^r 




9 f [{P - (l)a 9/3^ + {q- k)p + {k - p)^ g^/B] ^ 

_-g2 ^j^abej^cde {2g^^gp^ - g^^gp^, - ga|3g^^u) 



Additional Feynman rules for the model (47) 



P 



p 



2 " "mn ; 



a 



C 



r 
c' 



Xt 



m 




rvb 



m 



c 



2 S°'^^rt , 



p 



- 91,^"^, 



ky y a 



9{k + p), r'%t, 



- igf^^ a 



abc 

nm ' 



n 
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rv6 



a c „ ,, a Cy 
yr 



face xbde i fade £bce\ 



^6 ci"^ d-^ 



r 



a c. 



s 



^ -b rf^^ 



/ 



/ 

/' 



'-'b d'"" ''b d^"" +r'r'i26rsStu-Srt6su-6rJts)]. 

In addition to this: 

a) each closed loop brings a factor (27r)~^, 

b) each fermion or ghost loop gives an extra minus sign, 

c) arrows on the Majorana spinor lines should be ignored in calculating the symmetry 
factors. 

Dirac matrices in 4 — 2£: dimensions 

We use the metric g^^i, = (1, —1,-1, ...), (7^^ = 4 — 2e . 

[1^, = '2gf,u, 7t,l^ = 4 - 2e, 7m7i.7m = (2^ - 2)7,,, 

7M7i^7p7M = ^9up - 2e7i^7p, 7M7i^7p7a7M = 2e lulp^a - '^lalplv, (59) 

[75,7m]+ = 0' 75 = -!' tr75 = 0, tr/ = 4, tr(7^7^) = 4^1^^, 

^"l^ilplulalp) = K9,iu9cl3 - 9pa9uP + 9pl39ua), tr (7^^ . . .7^3^^ J = . 



The a- and /3-matrices of the model (47) 

These real antisymmetric 4x4 matrices have an explicit representation in terms of 
the Pauli matrices: 



, 1 



(Ji \ 2 _ / -0-3 \ z _ [ icr2 



" = ' -1 j ' " = Ua j ' « V .a, 



, ^ . i(T2 \ o2 ^ ( 1 \ ^ ^ ( -l(X2 

' ^ ia2 r ' \ -1 r ' I y ■ 



Their relevant properties are 



tra" = tr/3" = tr (a*^/?*) = 0, tr(a"a*) = tr(/3"/3*) = -45"* 



(60) 



(61) 



16 



The supersymmetric regularization used in section 6 implies S^^ — 3 + e giving rise to 
the following relations: 

= ^^^^ = -3 - £, aVo;'' = (1 + = (1 + £)/3*, (62) 

whereas the standard dimensional regularization prescribes 

(^'■'■ = 3, a'a'' = P^'P'' = -3, a' a'a'' = a\ = /3* . (63) 

Properties of the Euler F-function 



r{z + i)^zr{z), r(i) = r(2) = i, r{N + i)^N\, 



r(l + x) = exp 



n=2 



n 



(64) 



where 7 is the Euler constant and ( the Riemann function. We note that 7 and ({2) 
do not occur in )CTZ'G, and consequently in the renormalization group functions. 



One-loop integration formulas 

We choose a volume of the unit sphere in 4 — 2s dimensions to be |^ . 

J dp {p^Y — ^ 

dp in^ r{a + p-2 + e)r{2 - a - e) 



I 



dq 



i7r2r(l - e)T{a + 13-2 + e)T{2 -a- e)T{2 - /3 - e) 



(p2)«+/3-2+sr(a)r(/3)r(4 2^) 



(65) 
(66) 

(67) 



dq qf, _ m"^ p^ T{1 - e)V{a + 13-2 + e)T{3 -a- e)T{2 - 13 - e) 



q2a(^p _ 



(p2)"+/3-2+£r(a)r(/3)r(5 -a-P-2e) 



(68) 



dq q^,qu 



q2a(^p _ qyp 



iTT^ r(l - e)T{a + f3 - 3 + e)T{3 -a- e)T{2 - 13 - e) 
(p2)«+/3-2+^r(a)r(/3)r(6 -a- [3 -2s) 

x[{a + P-3 + e){3-a- e)p^p, + ^ (2 - ^ - e)g^y] (69) 



dq q^^quq\ _ itt^ T{1 - e)T{a + [3 - 3 + e)T{4: - a - e)T{2 - /3 - e) 



q2a^p _ qy0 



(p2)«+/3-2+er(a)r(/3)r(7 -a-/3-2s) 



X [{a + /3 - 3 + e){4 - a - e)p^PyPx + -{2 - j3 - e)p^{p^gv\ + Pvg„\ + Pxg^ci^)] (70) 
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Two- loop integration formulas [17] 



,2 \ a+/3+7+cr+p-4+2e 



[tTT 



2^2 



dt dq 



V{a,P,^,a, p) . 



V(a,l,7,l,l) 



r3(l - £)r(-l + 2e)r(l -a- e)V{l - 7 - s)V{a + 7 - 2 + 2s) 
r(a)r(7)r(3-a-7-3£) 
r(3-a-7-3£) r(a + 7-l + e) r(a) 



r(2-a-7-£) r(a + 7-2 + 3£) V{a-l + 2e) 

r(7) r(2-a-2e) r(2-7-2e)' 



+ 



r(7-l + 2£) r(l-a) 



r(i-7) 



(71) 



y(Q;,^,l,l,p) = 



X 



X 



r3(l - e)T{2 -a- e)T{2 - f3 - e)T{2 - p - s) 

r(2 - 2s)r(a)r(/3)r(p) 

(-)"r(n + 2 - 2£)r(m + n + a + /3 + p- 2 + 2£) 
^ m!n!(n + 1 - £)r(4 - m- Q!-/3-p - 3e)T(m + n + 2-e) 

,n=0 ^ ' ^ I I y \ / 

1 1 



E 



(n + p)(m + n + Q; + p — l + £) (n + p)(m + n + /3 + p — l + £ 
1 1 

+ 7 ^ T + 



m + n + Q;)(m + n + Q; + p— 1 + e) (m + n + /3)(m + n + /3 + p— 1 + 

1 1 

+ 



+ 



(m + n + Q;)(n + 2 - p - 2£) (m + n + ^)(n + 2 - p - 2e) 



(72) 



Individual two-loop integrals 

Here we write down the relevant integrals y(a, /3, 7, cr, p) with all the arguments 
being positive integers, retaining the ^, - and 0(1) terms. 



y(i,i,i,i,i) = 6C(3) 

1/(2,1,1,1,1) 
1^(1,1,1,1,2) 
1/(2,2,1,1,1) 
1^(2,1,2,1,1) 

y(2,i,i,2,i) 

1/(3,1,1,1,1) 



1 1 1 

2^2 ~ 2e ^ 2 
1 1 



1 _ 5 

£ ~ 2 
11^ 



£^ e 
1 5 11 

4£2 + 8£ ^ 16 
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Pole peirts of the essentially three-loop integrals of the form 

(A;2)3s-i f dp dq dt Yip.q.f.k) 



(i7r2)3 J p2q2t2(^k-p)^{k-q)^{k-t)^{p-qy{p-t)^{q-t)^ 



/ xs 2 61 877 4 , 

r= (P-tf ^ -3?-i8?-To8in«'^) 

xfi,2 1 41 31 6,,^, 

(p-f)V ^ ^((3) 

2 4Q 4 

e 

(k-qrip-tr ^ 



(k - q)V 



{k - q)Yk 



12£3 24^2 48£ 
1 65 865 



4e3 24^2 48£ 



3e^ 

1 3 53 
1 17 199 



6^3 12^2 24e 
,,,6 2 1 49 531 

^^-^^^ =^ 8i^ + 48i^+96i 

6^3 2^2 Qs 
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